Magnetic susceptibility of diluted pyrochlore and SCGO antiferromagnets 
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We investigate the magnetic susceptibility of the classical Heisenberg antiferromagnet with nearest- 
neighbour interactions on the geometrically frustrated pyrochlore lattice, for a pure system and in the 
presence of dilution with nonmagnetic ions. Using the fact that the correlation length in this system 
for small dilution is always short, we obtain an approximate but accurate expression for the magnetic 
susceptibility at all temperatures. We extend this theory to the compound SrCrg-g^^Gaa+gajOig 
{SCGO) and provide an explanation of the phenomenological model recently proposed by Schiffer 
and Daruka (Phys. Rev. B 56, 13712 (1997)). 

PACS numbers; 75.10.Hk, 75.30. Cr, 75.10.Nr, 75.40. Mg 

The study of magnetic systems with competing inter- 
actions has uncovered a wide variety of different physical 
phenomena M, such as glass transitions and the exis- 
tence low-temperature disordered phases. In a class of 
problems, which has received a great deal of attention 
recently B, competition between nearest-neighbour an- 
tiferromagnetic exchange interactions arises due to the 
fact that, in a group of g > 3 interacting spins, it is 
not possible for each spin to be antialigned with all its 
q — 1 neighbours. Geometric frustration of this kind can 
give rise to macroscopic ground-state degeneracies, and 
as early as 1956, a type of lattice was identified |^ for 
which the ground-state degeneracy is particularly large. 
In this kind of lattice, the frustrated units (triangles for 
the kagome and tetrahedra for the pyrochlore lattice. 
Fig. |l|) are arranged to share sites (neighbouring units 
share one spin) instead of bonds, as is the case in the 
more familiar triangular and face-centred cubic lattices. 

Since the manifold of ground states does not provide 
an intrinsic energy scale, any perturbation to the simple 
nearest-neighbour exchange Hamiltonian has to be con- 
sidered strong and can potentially select different low- 
temperature physics from the vast range of possibilities 
provided by the macroscopic degeneracy. One specatcu- 
lar example is the recently experimentally discovered 
magnetic analogue of ice [^-|| . The recent surge in theo- 
retical research on these systems results in large part from 
experimental developments (see, e.g., Refs. [Q,^^). Cur- 
rently, a systematic study of compounds is under way, 
in which different rare earth and transition metal ions 
are placed on the pyrochlore lattice, each of which comes 
with its own peculiar properties (such as anisotropic or 
longer-range interactions), so that the space of possible 
Hamiltonians is mapped out increasingly well. 

In this Letter, we study the magnetic susceptibility, Xj 
of a classical pyrochlore Heisenberg antiferromagnet as a 
function of temperature, both for the pure system and in 
the presence of disorder. The magnetic susceptibility is a 
particularly interesting quantity, since for a large class of 
geometrically frustrated magnets, its inverse surprisingly 
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FIG. 1. (a) The 'pyrochlore lattice, a network of cor- 
ner-sharing tetrahedra. (b) The trilayer SCGO lattice, con- 
sisting of corner-sharing tetrahedra and triangles. The top 
and bottom kagome layers are denoted by solid and broken 
lines, respectively, and the intervening triangular lattice by 
circles. By removing the spins labelled by crosses, only the 
encircled site of the q= 1 unit thus generated is left occupied. 



stays close to the linear Curie- Weiss law down to tem- 
peratures much lower than the Curie- Weiss temperature, 
Qcw^ where mean- field theory predicts a transition [EJ. 
(Eventually, at Tp <^ Qcw, such magnets usually, but 
not always |^, freeze S, a feature which is absent from 
the simple classical pyrochlore Heisenberg model P-O.) 
However, to our knowledge, no analytical expression, ex- 
act or approximate, which is valid at all temperatures, 
has so far been proposed for pyrochlore or SCGO antifer- 
romagnets (Fig. P, although for the kagome lattice, there 
has been some work at T=0 |13| ] , for infinite-component 
spins |14| , and using a high-order high-temperature series 
expansion |15| . 

Although the presence of at least a small amount of 
disorder is inevitable, its effect on pyrochlore magnets 
has so far only been studied qualitatively Pjl^]; again 
there has been some work for the kagome lattice at T = 
]l6[ . Recently, Schiffer and Daruka |17| have made the 
intriguing observation, described in detail below, that the 
deviation from linear Curie- Weiss behaviour at low tem- 



peratures for strongly frustrated magnets appears almost 
universally as a downturn of x~^, and that this tendency 
becomes stronger as the dilution is increased. This effect 
still awaits a theoretical explanation. 

We present here a theory which provides an expression 
for the energy and susceptibility of a classical Hciscn- 
berg pyrochlore antifcrromagnet which is exact at zero 
temperature and asymptotically correct for large tem- 
peratures, being a good approximation in between. Our 
treatment of the effect of dilution with non-magnetic ions 
applies in the limit of low disorder and is shown by Monte 
Carlo (MC) simulations to give reliable results for dilu- 
tions as large as 20%. We explain the empirical findings 
of Ref. 113 for SCGO by identifying two relevant mech- 
anisms responsible for the downturn of x~^j 0^6 due to 
dilution, the other present even in a pure system. 

We start with the fundamental observation that the 
spin-spin correlations of a pure classical Heisenberg an- 
tifcrromagnet on the pyrochlore lattice are always short- 
ranged [p|-[l2| . A small amount of dilution does not affect 
this property and induces neither ordering nor glassiness 
p. This ceases to be the case, at the latest, when 1/4 
of all sites are vacant. At this dilution, the possibility 
of moving local clusters of spins at no cost in energy is 
lost as the ground-state degeneracy of the magnet ceases 
to be macroscopic [|2[. We make use of this property 
by treating the weakly diluted lattice as an arrangement 
of units of g < 4 spins, which - given the small correla- 
tion length characteristic of this regime - are treated as 
if they were decoupled. 

For this approximation, which we refer to as the single- 
unit approximation, to be useful, we have to know the 
properties of the individual units. These can be ob- 
tained exactly using a Hubbard-Stratonovich transfor- 
mation [Q . We define the Hamiltonian for a group of q 
spins Si in the presence of a magnetic field B, to be 
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where we have dropped a constant in the last equality. 
The sum on (i,j) runs over all pairs. J > is the anti- 
ferromagnetic bond strength, and Lg is the total spin of 
the unit. The g- factor and the Bohr magneton have been 
absorbed into the definition of B. The partition function 
in zero field, Z^, is calculated to be 



Zq= /{df7,}exp(-/3JLV2) 



where the {Jl,} are the directions of the {8^} and /? = 
1/T is the inverse temperature. The final integrals may 
be evaluated in terms of the error function Erf, and we 
obtain, up to an overall constant: 



=^ 0.06 




FIG. 2. Susceptibility per spin, x, versus temperature from 
the single-unit approximation (line) and MC simulations on 
10976 spins (diamonds). Inset: Zero-field energy per spin 
versus temperature, both in units of J. 



Z2iT) = T[l-e-Ty (2) 

ZsiT) = ri {3Erf(V27r) - Erf(3V27r)}, (3) 

Z^{T) = T^ hETf(^/2/T] - ETf(y/8/T] 

- V778^(l-e"*) (3 + 2e-*+e-*)} (4) 

The susceptibility per spin, Xgy can be obtained from 
these expressions via the fluctuation-dissipation theorem: 
Xq = ((M2) - {Mf)/{3Tq). Noting that for one unit, 
AP = L^ is proportional to the energy, we finally obtain 
(using (M) = 0), X, = 2T/{iqJ)4r ^^Z,{T). 

For the full system, the Hamiltonian is given by 
the sum over the Hamiltonians of all N units, H ~ 
'Y^^^i i?,j" (B/2); here, q^ is the number of spins in unit 
a, and the magnetic field is divided by two to avoid dou- 
ble counting. 

Let us first consider the clean system, x = Q, where 
all units have g = 4. In the inset to Fig. H, we plot 
the energy per spin versus temperature in zero field, as 
obtained from the single-unit approximation and from 
Monte Carlo simulations. This plot demonstrates the 
quality of the approximation. The theoretical curves in 
Fig. H are scaled overall to take into account the fact 
that in the full system, each spin belongs to two tetra- 
hedra and has six rather than three nearest neighbours. 
Fig. m also shows a plot of the susceptibility. At high tem- 
perature, the analytic approximation and the simulation 
data both follow the Curie- Weiss law and therefore agree. 
As the temperature is lowered, the analytical result is 
only approximately valid, since it neglects the emerg- 
ing longer-range correlations. Still, since these are never 



strong, the disagreement never exceeds 5%. Remarkably, 
at the lowest temperatures, the two meet again, and at 
T — 0, our theory gives the exact result x(0) = 1/(8 J). 

Next, we turn to the diluted system. At high T, the 
effect of dilution is to decrease the average number of 
bonds of each spin, thus giving a reduced Curie- Weiss 
temperature Qcwix) = (1 — x)QcwiO) @- The main 
effect of disorder is seen by considering the susceptibil- 
ity at T = 0. On diluting the pyrochlore lattice, the 
probability of a tetrahedron containing q spins is given 

by Py{x) = (f\ (1 - a;)%'*-«. The effect of a magnetic 
field on the total spin of a unit is the same for all g > 2. 
From Eq. |l|, we can read off that -ffg(B) is minimised by 
choosing Lg — B/(2J), which is possible for all magnetic 
fields of strength B < 2Jq. A difference in the magnetic 
susceptibility per spin arises only because the total mag- 
netisation is shared between a total of g < 4 spins in a 
diluted unit. In any case, the susceptibility is finite since 
Lq oc B. This is in striking contrast to the case oi q — 1, 
where at T = an infitcsimal field suffices to align the 
spin, and hence Xi(0) is infinite. 

In the full system at low dilution, a unit with q = 1 
in general corresponds not to an isolated spin but to one 
which also belongs to a unit with q > 2 (see Fig. m. 
In this regime, there remain sufficiently many undi- 
luted tetrahedra to generate a macroscopically degener- 
ate ground state. This in turn provides the possibility of 
reorienting local spin clusters at no cost in energy p2[ . 
As a result, the correlation length remains short at all 
temperatures, and the spins in the q — I units, which are 
well-separated at low dilution, behave like paramagnetic 
spins even though they are not isolated. 

The behaviour at low temperatures is most easily dis- 
cussed in terms of the thermally averaged square of 
the total magnetisation of the system, ^M^(a;,T)). By 
equipartition, a unit with q > 3 has L^ « 3T/ J, and one 
with g = 2 has L2 « 2T/J, which vanishes as T — > 0. By 
contrast, a unit with q = I has Lf = S^ = 1. Since the 
magnetisation of the system equals the sum of the mag- 
netisations of the units, we obtain M = X]q=i L'-'^V^, 
where the factor 1/2 accounts for the fact that each 
spin belongs to two units. Therefore, (A/^(x,T)) — 

ELi (L^"^ ) V4 + Ea^c' (l^"^ • L(")' ) /4. Wc find that 

Ea^a' \ ■'-''"'' ' ■'-'''"■' ) = at low temperature, as would 
be appropriate for completely uncorrelated L^"^ . How- 
ever, our simulations suggest that this sum vanishes not 
because the individual terms are zero, but rather because 
the terms in the sum exactly cancel for T ^ 0. Hence, 
(M2(x,T)) = 6TNa~x)x{x,T) = E, 7V,(a;)L2(T)/4, 
where Nq = NPq{x) is the number of units with q spins, 
and L^ is given by the single-unit expression derived 
above. The estimate (M2(0)) = 7Vi/4 = Nx'^{l - x) 
is asymptotically exact in the limit a; ^ 0. For increas- 
ing dilution, terms higher order in x change the situation. 
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FIG. 3. A configuration containing two q = 1 units but 
witli vanishing magnetisation at low T. 



For example, the configuration depicted in Fig. g|, which 
is improbable at small x, contains two q = 1 units but 
has a vanishing magnetisation at T = 0. 

In Fig. |j, we plot the inverse susceptibility xi'T)^^ for 
different dilutions in the range < a^ < 0.2 at low tem- 
perature. Note that for x as large as 0.2, the agreement 
of our theory with the MC simulations is excellent. 

It is worth emphasizing that at any nonzero dilution, 
the low temperature susceptibility is dominated by the 
q — I units. This is analogous to a Griffiths-McCoy ef- 
fect as it arises from the rare, local event of three neigh- 
bouring sites being vacant. The temperature T below 
which the paramagnetic regions dominate can be defined 
as the point where the magnetisation due to the q — 1 
units equals that of the other units combined, which gives 
T (x x^ . This is why the temperature at which the down- 
turn of x^^ becomes visible is small and increases with 
disorder. 

Finally, we address in detail the work by Schiffer and 
Daruka jl^l, who proposed and successfully used a two- 
population model to fit the measured x to a form x — 

7=-^i l-7frr^ — . Here Ci and 0,„i are the Curie constant 

and Curie- Weiss temperature of a 'correlated' population 
which forms momentless clusters as T ^ 0, while C2 
and Qw2 are the parameters for an 'orphan' population, 
which was surmised to be excluded from the correlated 
population. Analysing experimental data on SCGO from 
Refs. [|9[, they found that 9^,2 can be set to zero, so 
that the orphan population in this case appears truly 
paramagnetic. For this compound, they present a series 
of different dilutions, independently determined to be in 
the range 0.11 < x < 0.61, and find that Qwi{x) = 
(1 — a;)0mi(O) and C2(x) ex x. 

The model we have presented above for the pyrochlores 
also applies to SCGO |12] , which consists of two kagome 
layers connected by an intervening triangular layer, and 
can be thought of as a slab of pyrochlore cut in a (111)- 
dircction (Fig. m. Every other kagome triangle is associ- 
ated with a site in the triangular layer, and these together 
form q — A units which generate the short correlation 
length required for our theory [ |l2| . 

The two main differences for a theory of SCGO are 
as follows. Firstly, only two spins (in a triangle not as- 
sociated with a site of the triangular lattice, as shown 
in Fig. ya) have to be removed to generate a. q = 1 
unit. Hence these units occur with a probability 3a;^/2 
at low dilution rather than Ax^ as in the pyrochlore. Sec- 
ondly, an additional mechanism is present even for an 
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FIG. 4. The inverse susceptibility for different dilutions x. 
The lines are given by the single-unit approximation without 
any adjustable parameters. Data points are from MC simu- 
lations of a lattice of 5488 tetrahedra. 



undiluted system, which generates a paramagnetic re- 
sponse in SCGO. If the interlayer coupling, which acts 
between spins of the triangular and the kagome layers, 
has strength J' different from J, the kagome intralayer 
coupling, the Hamihonian for a g = 4 unit reads -ff4(B) — 

(J/2) {(L' - B/{2J)f - (1 - J'/J)B -St}, where St is 
the spin in the triangular layer, and L' = (J'/J)Sf -I- 
^ j^ Si. In the zero-field ground state, each q = A unit 
has a finite magnetisation of magnitude 1 — J' /J. These 
individual magnetisation vectors are aligned at T = 
even by an infinitesimal field, yielding an infinite suscep- 
tibility. Arguments within the single-unit approximation 
along the lines of those given above predict a susceptibil- 
ity that diverges as 1/T at low temperatures. The value 
of J' /J is not known accurately, and there is no reason 
for it to be exactly 1. For a system known to be clean, the 
divergence of the susceptibility could be used to estimate 
this quantity. 

As a result, at zero and at low dilutions, the model of 



Ref. 1 17 1 is reasonable even though the orphan popula- 
tion is not excluded from the correlated one but present 
even for the pure system and enhanced by the creation 
oi q = 1 units at small dilution. (Genuinely isolated 
spins of course appear at high dilution.) The empirical 
fit C2 (x) oc 2; is clearly at variance with our theory, since 
it yields C2(x) = (f - J'/J)^/21 + x^/U for smah x. 
However, there is no contradiction, since the experimen- 
tal data consists of only five points, and there are no data 
points sufficiently close to a; = to distinguish between 
the two functional forms of C2 . 

A quantitative comparison of our theory with experi- 
ment may be possible in the pyrochlores, since some com- 
pounds can be grown with very small amounts of disor- 



der 1 20 1 . Introducing controlled amounts of dilution with 
nonmagnetic ions can then produce a series of samples 
allowing measurement of susceptibility as a function of 
disorder and temperature. 
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